Migration of nuclear shell gaps studied in the d(^''Ne,p7)^^Ne reaction 
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The transfer of neutrons onto '^''Ne has been measured using a reaccelerated radioactive beam of 
Ne to study the (d,p) reaction in inverse kinematics. The unusual raising of the first 3/2+ level m 
■^*Ne and its significance in terms of the migration of the neutron magic number from N=20 to N=16 
is put on a firm footing by confirmation of this state's identity. The raised 3/2+ level is observed 
simultaneously with the intruder negative parity 7/2~ and 3/2~ levels, providing evidence for the 
reduction in the N=20 gap. The coincident gamma-ray decays allowed the assignment of spins as 
well as the transferred orbital angular momentum. The excitation energy of the 3/2+ state shows 
that the established USD shell model breaks down well within the sd model space and requires a 
revised treatment of the proton-neutron monopole interaction. 

PACS numbers: 21.10.Hw21.10. Jx23.20.Lv25.60. Jc27.30.+t29.38.Gj 



The monopole part of the nucleon-nucleon interaction 
is now recognised as having a major effect on nuclear 
shell structure far from stability [l[ Q . The interaction 
between valence protons and neutrons is sufficient to al- 
ter the energies of single-particle levels so that different 
magic numbers (or shell gaps) appear and this can sub- 
stantially affect the collective [§] and magnetic Q proper- 
ties and basic quantities such as the lifetime . Nucleon 
transfer reactions induced by light ions are an established 
experimental tool for studying single-particle structure 
0. Here we employ the (d,p) reaction in inverse kine- 
matics to explore the disappearance of the N=20 magic 
number (and its replacement by N=16) in the neutron- 
rich neon isotones. As will be shown, the measurement of 
the differential cross sections of the light ejectiles plus the 
coincident gamma decays of the residual nucleus brings 
a new power to this type of study. 

Recent work using other techniques has provided ev- 
idence for the emergence of N=16 as a magic number 
in this region, but has not identified the single-particle 
structure in an unambiguous manner through measure- 
ments of the spectroscopic factors and spins. In a study 
of the /3-decay of ^^F Q the increased energy of the 
neutron orbital was inferred. This made use of a pre- 
liminary analysis of the present work Q , and concluded 
that the energy shift was consistent with the monopole 
effect 0. In a study of ^''Ne using the (d,p7) reaction 
but without detecting the protons [91] , a reduced gap be- 
tween the 0c?3/2 and higher negative parity orbitals was 



deduced. This agreed with nucleon removal studies [lOj. 
Finally, in recent studies of ^'^O by transfer llj and ^^O 
by proton removal [13] the 0^3/2 state was found to have 
an increased excitation energy but the required modifi- 
cations to the shell model interaction were not mutually 
consistent [ff|, [f^] . While an extensive review including 
the emergence of the N=f6 magic number has recently 
been published |2], further quantitative data are needed 
in order to understand this monopole effect properly. 

In this spirit, we have investigated the levels of ^^Ne, 
in which the valence neutron can occupy the fsi/2 or 
0^3/2 orbitals, or the higher lying negative parity orbitals 
(0/7/2, f]53/2, • ■ •)• This explores the N=16, N=20 and 
N=28 gaps. Existing information about levels in ^^Ne 
is included in Table I. The /3-decay work showed 0], in 
addition to the two excited states listed, three other pos- 
itive parity states above 3.3 MeV including a state seen 
also in nucleon removal from ^^Ne and proposed to have 
,r = 5/2+ [13]. The present 3.33 MeV state was pro- 
posed [2J] to have negative parity (along with the 4.07 
MeV state) based on the ("C,i^O) selectivity. 

A pure beam of ^^Ne ions, at 10.6 -A MeV and 2 x 10^ 
particles/sec, was provided by the SPIRAL facility at 
GANIL. The beam impinged on a self-supporting target 
of deuterated polythene (CD)2 with thickness I mg/cm^. 
The target was surrounded by TIARA fSj, [ij] , an array 
of silicon detectors spanning 85% of 47r. The TIARA cov- 
erage from 36° to 144° comprised an octagonal "barrel" 
of position sensitive (resistive strip) detectors, 400/im 
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thick. The most backward angles (144° to 169°) were 
covered by an array of 400/L(ni thick double-sided Si strip 
detectors (DSSSDs). In very close proximity to the tar- 
get (50mm) the front faces of four EXOGAM segmented 
Ge detectors [31 formed four sides of a cube spanning 
66% of 47r. The full-energy peak efficiency for 7-ray de- 
tection was estimated to be 15% at 1.0 MeV. During the 
experiment, however, the Ge detectors were in operation 
for only a fraction of the beam exposure. Beam particles 
and beam-like reaction products emerging near zero de- 
grees were analysed by the VAMOS magnetic spectrom- 
eter [1^. The focal plane detectors gave Z identification 
by AE-E and mass by time-of- flight [l3]- The start sig- 
nal for the timing was provided by a thin (lO/itm) plastic 
scintillator mounted 0.5m upstream of the reaction tar- 
get. The data acquisition was triggered if an event was 
recorded in any of the TIARA silicon detectors. Any co- 
incident event in VAMOS was recorded and, in that case, 
also any further coincident event in EXOGAM. Only the 
core signals p^] from EXOGAM were available, limiting 
the accuracy of the Doppler energy corrections for 7-rays 
emitted by the beam-like particles {(3 w 0.1) to a FWHM 
of 65 keV at 1 MeV. 

Products of ^"^Ne reacting with ^^C in the target were 
evident in the TIARA singles spectra of energy against 
angle, comprising knock-on ^^C plus protons and a- 
particles from fusion-evaporation, but were eliminated by 
requiring a beam or beam-like particle in VAMOS. The 
elastically scattered deuterons were detected just forward 
of 90° [17| in coincidence with scattered beam recorded 
in VAMOS. These measurements provided an absolute 
normalization of all the differential cross sections. Pro- 
tons from the (d,p) reaction were recorded from the most 
backward angles until the increasing energies resulted in 
their penetrating the barrel near 95°. The protons corre- 
sponding to the most forward cm. angles were detected 
in the DSSSDs between 135° and 170°, with improved 
energy resolution and lower energy thresholds compared 
to the resistive strip detectors in the barrel. 

The particles observed in TIARA in coincidence with 
^^Ne were used to calculate the ^^Ne excitation energy, 
according to the measured energies and angles and the 
kinematics for (d,p). The resolution in excitation en- 
ergy was strongly angle dependent, owing to the target 
effects and the detector resolution 18|. In total, 8 an- 
gular bins were used, including 5 in the barrel. Fig. [T] 
shows example excitation energy spectra. These were fit- 
ted with a series of Gaussian peaks, using the 7-ray data 
to identify the peaks to be included and their energies. 
The 7-ray energy spectrum from all triple coincidence 
^^Ne-p-7 events, after Doppler correction, is displayed in 
Fig. [2ja). Gamma-ray peaks were observed at energies 
of 1.68, 2.03, 2.35 and 3.33 MeV (± 0.04 MeV). As dis- 
cussed below, the 7-ray at 2.35 MeV is in cascade with 
that at 1.68 MeV, suggesting a state at 4.03 MeV. With 
these 4 excitation energies and the ground state fixed, 
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FIG. 1: Excitation energy spectra deduced from protons de- 
tected at laboratory angles of: (a) 109°, (b) 125°, (c) 149°, 
(d) 164° (color online). 



the spectra for the various angle bins were each fitted 
with 5 peaks of equal (but angle-dependent) width. 
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FIG. 2: Doppler-corrected gamma-ray energy spectra 
recorded in coincidence with (a) all protons and (b) protons 
populating the state in ^'^Ne at 4.03 MeV; (c) levels in ^^Ne. 

The angular distributions from the (d,p) reaction pop- 
ulating states in ^^Ne are shown in Fig. [3l together with 
calculations using the Adiabatic Distorted Wave Approx- 
imation (ADWA) with standard parameters that 
have been shown to work well for comparisons with large- 
basis shell-model calculations [^. The nucleon- nucleus 
potential of Bechetti-Greenlees was employed The 
elastic (d,d) cross section was measured [l3| between 
25° and 45° (cm.) simultaneously with the (d,p) and 
fitted using optical potential parameters obtained from 
d+26Mg at 6 -A MeV [2lj. This calibrated the product 
of beam exposure and target thickness. 

The angular distribution for transfer to the ^^Ne 
ground state exhibits a characteristic dip near 135° (17° 
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in the cm. frame) indicating an ^ = transfer. The an- 
gular distributions of the 1.68 MeV and 2.03 MeV states 
are consistent only with 1 — 2 transfer and are clearly 
different from that of the ground state. The other strong 
peaks at 3.33 and 4.03 MeV are in a region where USD 
shell-model calculations (such as those included in Ta- 
ble I) predict little transfer strength to positive parity 
states. For the 3.33 MeV state, the fall-off in yield with 
increasing cm. angle (decreasing laboratory angle) is too 
rapid for high ^-transfers and implies i ~\. For the 4.03 
MeV state, the experiment provided only a lower limit 
on the cross section at some angles. This resulted from 
the loss of some counts below the energy threshold of the 
barrel detectors, which is increasingly important near the 
end of the barrel (135°) owing to the kinematics. The 
observed yield near 90° is sufficient, however, to rule out 
any rapid fall-off with increasing cm. angle, and in fact 
the data are consistent only with an ^-transfer of at least 
^ = 3, with i — 2 giving too low a yield near 90° com- 
pared to that near 170°. Since there are no £ = 4 orbitals 
nearby in energy in the shell model, an £ = 3 transfer is 
deduced for the 4.03 MeV state. 

The deduced ^-transfers can be extended to spin as- 
signments using the observed 7-ray decay scheme. The 
7-ray spectrum obtained using a restricted gate on exci- 
tation energy near 4 MeV is shown in Fig. El^b). When 
the quadruple-coincidence ^^Ne-p-7-7 data are analyzed, 
the limited statistics are sufficient to demonstrate that a 
gate on the 1.68 MeV peak highlights the peak at 2.35 
MeV and, similarly, a 2.35 MeV gate highlights the 1.68 
MeV transition Il4[. Such a cascade is of course expected 
for a (5/2^, 7/2^) £ ~ 3 state, decaying via an £ = 2 
state, and the cascade from a 7/2~ state would select just 
the 5/2"*" state. From simple shell model considerations, 
the lowest £ — 3 state is expected to be 7/2^. Indeed, 
the isotonic reaction ^^Mg(d,p)^''Mg [2l!| strongly popu- 
lates an £ = 3 state at 3.76 MeVjwhich is assigned 7/2^ 
and decays via the 5/2+ state |22|. The combined evi- 
dence implies that the 4.03 MeV state in ^^Ne should be 
assigned 7/2~ and also that the intermediate 1.68 MeV 
state in the decay is 5/2+. Furthermore, the reaction 
yield to the 1.68 MeV state is lower than that for the 
2.03 MeV state in ^^Ne. This reinforces the 5/2+ assign- 
ment, since the 5/2+ hole state (naively i^(ls^^2 ® 0*^5/2 )) 
should be populated more weakly than the 3/2+ parti- 
cle state (1^(0^3/2)), as also observed in ^''Mg By 
elimination, the £ = 2 state at 2.03 MeV can be inferred 
to be the 3/2+ state expected from systematics and shell 
model considerations. These arguments for the £ = 2 
states are all supported by high-energy nucleon removal 
[13] and /3-decay studies. The 9Be(26Ne,25Ne) reac- 
tion was seen to populate the 5/2+ hole state but not the 
3/2+ state The ^^F ;3-decay was found to favour the 
2090 keV state relative to that at 1702 keV and compari- 
son with USD calculations suggested that the lower state 



must be 5/2+, leaving the higher energy state as 3/2+ 
01 . Finally, the lowest energy 1=1 state is expected 
to be 3/2~ in the shell model, and hence the 3.33 MeV 
state in ^^Ne will correspond to the strongly populated 
3.56 MeV 3/2~ state in ^''Mg. As expected, this state in 
^^Ne 7-decays to the 1/2+ ground state. 

The excitation energies and deduced spin-parity as- 
signments are listed in Table I. The adopted energies for 
the first two excited states are defined to 1 keV precision 
by the beta-decay results 0] whilst the weighted average 
energies for the negative parity states are 3325 ± 25 keV 
and 4055 ±25 keV. The precision in the present work is 40 
keV owing to the limited statistics and the large widths 
of the gamma-ray peaks. The relative ordering of the 
lowest 7/2~ and 3/2 ~ levels is the same as in the isotone 
^''Mg, but they are further apart in ^^Ne. This contra- 
dicts Monte-Carlo SDPF-M shell model predictions for 
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FIG. 3: Differential cross sections for final states in ^^Ne, 
labefed by excitation energy and compared with normalized 
ADWA calculations. The solid line for each state corresponds 
to the adopted ^-transfer and spin as included in Table I. 
For the points indicating lower limits (arrows), the detector 
energy thresholds reduced the efficiency. 
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TABLE I: Results for states in ^Ne. The measured excitation energies are compared with previous work, where the experimental 
errors are 40 keV (present), of order 1 keV 0, [ISli 30 keV [23l and 50-80 keV [25l. States observed in the beta-decay work but 
not expected in transfer have been omitted, as have states at 4.7±0.1 MeV j25l] and 6.28±0.05 MeV [24| . 
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^^Ne [101 ■ Much more dramatically, whilst the 5/2+ state 
is at a similar energy in the two nuclei (reflecting in part 
the separation of the 0(^5/2 and lsi/2 orbitals), the low- 
est 3/2+ state moves from about 700 keV below the 5/2+ 
in ^"^Mg to 250 keV above in ^^Ne. This is unexpected 
according to standard USD shell-model calculations (26| 
(Table I) and reflects a raising in energy of the 0^3/2 or- 
bital in neutron-rich nuclei. This can be reproduced with 
a revised treatment of the monopole interaction 0, [l^l . 
The present data show that the raising of the 0^3/2 or- 
bital simultaneously opens the N=16 gap and closes the 
N=20 gap as the negative parity states are at similar 
energies in ^^Mg and ^^Ne. 

Having adopted the assignments listed in Table I, 
the ADWA calculations can be scaled to the experimen- 
tal data and spectroscopic factors (5*) deduced. These 
are also listed in Table I, where a value of unity would 
indicate transfer to an unblocked pure single-particle or- 
bital. The uncertainty arising from the theoretical model 
is estimated as 20% %. The 1/2+, 3/2" and 7/2" states 
all carry a substantial part of the single-particle strength. 
The 5 /2+ spectroscopic factor is smaller because it is pre- 
dominantly a hole state and as such effectively measures 
the amount of 1^(0^5/2)^ to v{lsi/2)^ pair excitation in 
^^Ne. The 3/2+ state is evidently of complex structure; 
besides the 0^3/2 single-particle structure it includes, for 
example, a component of ^'*Ne(2+) ® lsi/2- According 
to USD calculations, the remaining 0^3/2 single-particle 
strength is distributed widely, with much in the third 
and fifth 3/2+ states between 4.5 and 6 MeV. We note 
that in the N=15 isotone ^'^O, the 3/2+ state was ob- 
served [llj (with S = 0.5±0.1) at an excitation energy 
well above that predicted by the USD interaction [26| . 
Interestingly, while this could be reproduced using the 
revised USD interactions [13], they fail to explain the 
excitation energy of the corresponding state in ^^O 12 1. 



In conclusion, the spins of the low-lying levels in ^^Ne 
have been firmly assigned. The inversion of the 3/2+ and 
5/2+ levels represents a significant breakdown of the USD 
model, which can be explained in terms of a monopole 
shift in the effective 0^3/2 single-particle energy 0, [l^] 



so as to open a shell gap at N=16. The measured spec- 
troscopic factors for ^^Ne provide the means to test in 
detail the shell model modifications invoked to explain 
this shift. The present data also indicate that the open- 
ing of the N=16 gap is accompanied by a closing of the 
N=20 gap . This confirms the inference from ^^Ne stud- 
ies |l^. Finally, the experimental methods of (d,p7) 
and (d,p77) employed here are likely to find wide ap- 
plicability in the exploration of nuclear structure with 
radioactive beams. 
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